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The Durham/UKST Galaxy Redshift Survey - VII. 
Redshift-space distortions in the power spectrum 
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ABSTRACT 

We investigate the effect of redshift-space distortions in the power spectrum parallel 
and perpendicular to the observer's hue of sight, P'^(fcj[,kj^), using the optically se- 
lected Durham/UKST Galaxy Redshift Survey. On small, non-linear scales anisotropy 
in the power-spectrum is dominated by the galaxy velocity dispersion; the 'Finger of 
God' effect. On larger, linear scales coherent peculiar velocities due to the infall of 
galaxies into overdense regions are the main cause of anisotropy. According to gravi- 
tational instability theory these distortions depend only on the density and bias pa- 
rameters via /3 K, Vl^ /h. Geometrical distortions also occur if the wrong cosmology 
is assumed, although these would be relatively small given the low redshift of the 
survey. To quantify these effects, we assume the real-space power spectrum of the 
APM Galaxy Survey, and fit a simple model for the redshift-space and geometrical 
distortions. Assuming a flat f7,„ = 1 universe, we find values for the one-dimensional 
pairwise velocity dispersion of Cp = 410 ± 170 km s""'^, and (5 = 0.38 ± 0.17. An open 
Vim — 0.3, and a flat ^m — 0.3, f^A = 0.7 universe yield Cp — 420 km s^^, (3 — 0.40, 
and (jp — 440 km s^^, j3 = 0.45 respectively, with comparable errors. These results are 
consistent with estimates using the two-point galaxy correlation function, ^(ct, tt), and 
favour either a low-density universe with fim ~ 0.3 if galaxies trace the underlying 
mass distribution, or a bias factor of 6 ^ 2.5 if f7m = 1. 
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1 INTRODUCTION 

When mapping the three-dimensional distribution of galax- 
ies redshift is generally the easiest, and often the only dis- 
tance estimator available. Redshifts, however, do not reflect 
real galaxy distances because of peculiar velocities relative 
to the Hubble flow. Whilst standard cosmologies assume 
that the Universe, and hence galaxy clustering, is isotropic 
in real-space, the effect of peculiar velocities and incorrect 
geometric assumptions can distort the picture in redshift- 
space. The anisotropy is dominated by the galaxy velocity 
dispersion in virialized clusters on the smallest scales. This 
is the 'Finger of God' effect seen in redshift surveys. The 
main cause of anisotropy on larger, linear scales is the co- 
herent peculiar velocity fleld of large-scale structure. The 
infall of galaxies into overdense regions makes their distri- 
bution appear slightly squashed in redshift-space, increasing 
the apparent overdensity. According to linear gravitational 
instability theory these distortions depend only on the den- 
sity and bias parameters via the combination /3 ~ Q.^ /h 
(Kaiser 1987). Incorrect geometric assumptions through the 
assumed cosmological model also add an anisotropy on all 
scales. These redshift-space distortions can be exploited in 



order to constrain properties of the galaxy velocity fleld, and 
fundamental cosmological parameters such as Q, and A. 

Redshift-space distortions can be investigated using the 
spatial two-point correlation function or its Fourier trans- 
form, the power spectrum. In an ideal situation these two 
approaches are equivalent. When applied to finite, noisy 
data, however, they have different strengths, and should be 
viewed as complementary. In particular the power spectrum 
approach has the advantage that uncertainty in the mean 
galaxy density only affects the fe = mode, whereas uncer- 
tainty in this quantity introduces errors in the correlation 
function measurement at all scales. Also, the power spec- 
trum is more naturally related to theoretical model predic- 
tions, and the modelling of these anisotropics is simpler in 
redshift-space. 

Early attempts to investigate the anisotropics in the 
power spectrum (Cole, Fisher, & Weinberg 1994) focussed 
on the linear regime, developing a simple estimator for j3 
by considering the quadrupole and monopole angular mo- 
ments of the power spectrum. This model was improved to 
account for a small-scale velocity dispersion (Cole, Fisher, 
& Weinberg 1995), but still ignored geometric distortions 
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due to cosmological assumptions. Any geometric distortion 
has a similar effect as that of /3, and if ignored this would 
lead to systematic errors. Whilst working at very low red- 
shifts, where the geometric distortions are negligible, this 
approach is valid, but at the higher redshifts probed by the 
future large galaxy and quasar redshift surveys, cosmolog- 
ical effects need to be considered. Reducing the full two- 
dimensional information available in P {k^\,'k±) to just the 
first two angular moments would make it impossible to dis- 
criminate the subtle angular differences between geometrical 
and redshift-space distortions. Ballinger, Peacock, & Heav- 
ens (1996) developed a model for P (A:||,kx) taking all of 
these effects into account. 

In this paper we investigate the redshift-space distor- 
tions in the power spectrum parallel and perpendicular to 
the observer's line of sight, P (fc||,kx), using the optically 
selected Durham/UKST Galaxy Redshift Survey. We apply 
the model from Ballinger et al. (1996) to the power spectrum 
obtained and derive estimates for f3, and the one-dimensional 
pairwise velocity dispersion, ap. Finally, we compare our re- 
sults to those obtained from the same survey through anal- 
ysis of the two-point correlation function (Ratcliffe et al. 
1998a). 



2 POWER SPECTRUM ANALYSIS 

2.1 Sample definition 

The power spectrum analysis is carried out on the volume 
limited sample, with z < Zmax = 0.06, of the optically se- 
lected Durham/UKST Galaxy Redshift Survey, described 
by Hoyle et al. (1999). Full details of the construction of the 
survey can be found in Ratcliffe et al. (1998b). 

2.2 Power spectrum estimation 

The power spectrum is estimated using the method out- 
fined in Hoyle et al. (1999) and Tadros & Efstathiou (1996). 
The Durham/UKST Survey is divided into four regions in 
order to apply the distant observer approximation (see sec- 
tion 2.c). Each region is embedded into a larger cubical vol- 
ume, which is rotated such that the central line of sight lies 
along the same axis of the cube in each case. The density 
field is binned onto a 256^ mesh, using nearest gridpoint as- 
signment. The power spectrum of each region is estimated 
using a Fast Fourier Transform, and the average of the re- 
sulting power spectra is taken. The result, binned logarith- 
mically into fc|| and kx is plotted in figure hi 



2.3 Survey geometry 

The Durham/UKST Survey covers a region of approxi- 
mately 20° X 75°. In order to retain information about 
the line of sight, to calculate P (fc||,kx), the data have to 
be divided into subsamples that subtend a small solid an- 
gle on the sky, and the distant-observer approximation ap- 
plied. The data are split into 4 regions each approximately 
20° X 20° during this analysis, and the average of the re- 
sulting power spectra is taken. The measured P (fcy , kx) is 
in fact a convolution between the true power spectrum and 
a function related to the sample window introduced due to 




Figure 1. P^(fc||,kx) estimated from the Durham/UKST sur- 
vey. Filled contours of constant logP are shown as a function of 
fcy and kx- Overlaid is the best-fit model for P^(fc|| , kx), assum- 
ing a flat flrn = 1 universe, with /3 = 0.38 and crp = 410 km s~^, 
and the real-space APM power spectrum. 



the assumption of parallel lines of sight. This makes the 
measured power spectrum more isotropic, and hence causes 
a small, yet systematic underestimate of /3. Although it is 
possible to correct for this effect (see Cole et al. 1994), the 
error in this approximation for a 20° x 20° window is < 1% 
and considerably smaller than the statistical errors obtained 
in this analysis. Hence for the purposes of this analysis it is 
ignored. 

The measured power spectrum is also convolved with 
the power spectrum of the window function. As we are con- 
sidering a volume limited sample, each galaxy carries equal 
weight, and hence the survey window function, VK(x), simply 
takes a value of unity in the volume of the universe included 
in the survey, and zero elsewhere. This is approximated using 
a catalogue containing a large number of unclustered points 
with the same radial and angular distribution as the survey, 
and its power spectrum is calculated in a similar manner to 
that of the data. Figure ti shows the power spectrum of the 
window function, IW"!^, as a function of fcy and kx. 

We can probe larger scales along the line of sight than 
across it; the depth of the survey is ~ 175/i~^Mpc, cor- 
responding to fc||mi„ ~ 0.04 h Mpc~^, whereas the maxi- 
mum scale probed at this depth across the line of sight in a 
20° X 20° region is ~ 85/i~^Mpc, corresponding to |kxmin| 
~ 0.07 h Mpc~^. The effect of this can be seen in the 
anisotropy of | PFp on the largest scales. Power on the largest 
scales across the line of sight is lost when the data are split 
into sub-samples. The choice of 4 sub-samples is a compro- 
mise between minimizing the effects of the distant- observer 
approximation and maximising the volume of k-space we 
can probe. 

The window function is sharply peaked in k-space, vary- 
ing as a steep power law, ~ k~^ or steeper; much steeper 
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Figure 2. The power spectrum of the window function. Filled 
contours of constant log{|l4^(A;||, kx)P) are shown as a function 
of fc|| and kx- For wavenumbers with fcy > 0.06 h Mpc~^ or |kx| 
> 0.1 h Mpc~^, the window function power spectrum is a steep 



power law, varying as ■ 



or steeper. 



than the expected galaxy power spectrum. Therefore the 
convolution should have little effect at wavenumbers where 
\W\^ takes a value less than about a tenth of its maximum 
value. At wavenumber separations where \W\^ takes a value 
greater than this, the convolution will have an effect on our 
estimate of the galaxy power spectrum. To avoid this, we 
consider only those regions with fcy > 0.06 h Mpc~^ or |kx 
> 0.1 /i Mpc~\ 



2.4 Mock catalogues 

In order to test our power spectrum estimator, and to esti- 
mate the statistical errors on the recovered power spectrum, 
mock catalogues are used. A total of 40 mock catalogues 
are constructed using the Hubble Volume rCDM simulation 
from the Virgo Consortium (Colberg et al. 1998). They have 
the same radial and angular selection as the Durham/UKST 
Survey, and approximately the same clustering amplitude. 
For further details on the construction of the mock cata- 
logues see Hoyle et al. (1999). The power spectrum of each 
catalogue was measured in exactly the same way as for the 
data. 

The fractional errors, calculated from variance in the 
power spectrum measurements across the 40 mock cata- 
logues, are displayed in figure H. The errors on the power 
spectrum estimate are lowest at high values of |kx|. This is 
because kx represents two dimensions in Fourier space, and 
hence the Fourier space volume sampled in any bin oc Afcy 
x|kx| A|kx|. 

As well as using the mock catalogues to estimate the 
statistical errors, we can also investigate possible systemat- 
ics, such as the effect of the window function, by fitting the 
redshift-space distortions in the mock power spectra, using 



Figure 3. The fractional errors in P^ {ku, kx), as estimated using 
mock catalogues, arc shown as a function of fcii and kx- 



the method described in the next section. The average of 
the 40 mock catalogue power spectra is shown in figure W, 
together with the best fitting model. The errors are con- 
siderably smaller than those of the data, as this is the av- 
erage over 40 realizations. The discrepancy on the largest 
scales is partly due to the effects of the window function. 
To minimise the effects of this uncertainty, the region with 
fc|| > 0.06 h Mpc"^ and jkx| > 0.1 h Mpc"^ is excluded 
from the fit. Also, fractional errors decrease with increasing 
kx I , as shown in figure H, and so the agreement between the 
mock catalogues and the model is best at large |kx|. The 
best-fit model for P^^^,^{k\\,'kj_), assuming a flat Qm = 1 
universe, has /3 = 0.37 and ap = 540 km s"'^. The pair- 
wise velocity measurement is in remarkable agreement with 
that found in the simulation; Jenkins et al. (1998) measured 
ap ~ 530 km s~^ for the whole simulation on large scales. 
The value of /3 depends on the precise bias prescription used, 
but is in good agreement with other estimates. Hence we 
conclude that our method produces a robust measurement 
of the redshift-space distortions, and any systematic uncer- 
tainties in our measurements, especially due to the effects of 
window function, are small on the scales considered. 



3 MODELLING REDSHIFT DISTORTIONS 

The model that we fit to the data is based on that presented 
in Ballinger, Peacock and Heavens (1996). For clarity and 
convenience we review the main details below. 



3.1 Geometrical distortions 

The power spectrum analysis was carried out assuming an 
EdS universe. If the true geometry is different then our dis- 
tance calculations will be wrong by a factor f± perpendic- 
ular to the line of sight, and /y along the line of sight (as 
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Figure 4. The average P^ . (fc||,kx) estimated from 40 mock 
catalogues, constructed from the Virgo Consortium's Hubble Vol- 
ume rCDM simulation. Filled contours of constant logP are 
shown as a function of fcii and kx- Overlaid is the best-fit 
kx), assuming a flat Vim = 1 universe, with 



model for P^ , (k 

mock ^ 



/3 = 0.37 and op = 540 km s"!. 



defined in Ballinger et al. 1996). Tlius we can define a geo- 
metric flattening factor: 



f± 4' 



+ 2nA)z + o{z^ 



(1) 



(BaUinger et al. 1996). For an open {Qm = 0.3, ^a = 0.0), 
and lambda-dominated (f2m = 0.3, f^A — 0.7) universe, F ~ 
1.01 and 1.03 respectively at the survey limit. As the deptii 
of the survey is increased this effect becomes much more 
important. The effect this has on the power spectrum is 
given by 



assume 



d(fc||,kj 






1 + M' 



(i^-O 



(2) 



k«/k, and n is the spectral 



(Ballinger et al. 1996) where /i 
index of the power spectrum. 

3.2 Redshift distortions 

On large, linear scales coherent peculiar velocities due to the 
infall of galaxies into overdense regions are the main cause of 
anisotropy. Kaiser (1987) showed that this anisotropy takes 
a very simple form in redshift-space, depending only on the 
density and bias parameters via the combination f3 ~ O!^ /b: 

P^(fc||,kx) = P^(fc)[l-f/3MT (3) 

where P^ , and P^ refer to the redshift-space and real-space 
power spectra respectively. The Kaiser formula assumes the 
distant-observer approximation, and hence is valid in thi s 
case due to the method of estimation of P (fc|| , k±) (§ 2.3). 
On small, non-linear scales anisotropy in the power- 
spectrum is dominated by the galaxy velocity dispersions in 



Table 1. Minimum x^ fits (over 23 bins) for Op and /3 to 
P (k\i , kx) estimated from the Durham/UKST survey. The open 
and lambda cosmological models both assume Q,m = 0.3, and 
SIa = 0.0 and 0.7 respectively. 



Cosmology o-p/km s 



P 



X 



EdS 


410±170 


0.38±0.17 


14.18 


Open 


420±150 


0.40±0.19 


14.11 


A 


440±200 


0.45±0.19 


13.78 



virialized clusters. This is modelled by introducing a damp- 
ing term. Assuming an exponential line of sight pairwise 
velocity (crp)n distribution (Davis & Peebles 1983), this is 
given by a Lorentzian factor in redshift-space: 



D [kiiOp] = 



1 -I- I (k/iap) 



(4) 



This model is only a good approximation in mildly non- 
linear situations. Cole et al. (1995) used numerical sim- 
ulations to show that it is effective for wavelengths A ^ 
15/i~^ Mpc, or fc ^ 0.4. At such wavenumbers assuming a 
Gaussian velocity dispersion has a negligible difference on 
the resulting model (Ballinger et al. 1996). 

Combining these effects leads to the final model: 



P^(fc||,ki 



P«(fc) 



1 I 2 

1 -f A* 



l + M 



i^-^) 



(3+1 

F2 



D Ikfiapl 



(5) 



(Ballinger et al. 1996) where a'p = cFp/ f\\ . There are several 
free parameters in the model; F, (3, Op, and P^{k), the under- 
lying real-space power spectrum. To reduce the uncertain- 
ties on the parameters of interest we assume the real-space 
power spectrum of the APM Survey galaxies (Baugh & Efs- 
tathiou 1993). When fitting to the rCDM mock catalogues, 
we assume instead a rCDM real-space power spectrum. 



3.3 x^ fitting 

We fit the model to the data using a x^ technique: 



X 



\ '^ \J^obs ^mod) 






(6) 



'p(k) 



The errors, crp(fc), are those estimated from the mock cat- 
alogue power spectra. Only those wavenumbers with fcy > 
0.06 or |kx| > 0.1, and k < 0.4 are used in the fit. The 
former constraints are applied to remove the effects of the 
window function, and the latter to prevent excessive non- 
linearity, where the model breaks down. 

An underlying cosmology is assumed and the l3 and 
cTp parameters are adjusted until a minimum value of x^ 
is obtained. The errors are calculated using Monte Carlo 
simulations of data sets drawn from the best fltting model, 
using the mock catalogue errors. The results are displayed 
in table m. 



* In power spectra crp is implicitly divided by Hq and quoted in 
units ft-l Mpc. (Ho = lOOh km s'^ Mpc'^) 
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density universe with Qrn ~ 0.3 if galaxies trace the under- 
lying mass distribution, or a bias factor of 6 ~ 2.5 if Qm ~ 1. 
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Figure 5. Contours of increasing x^ in the Op - /3 plane, assuming 
a flat Qrn = 1 universe. The solid contour represents the l-cr error 
in the Durham/UKST survey estimate, the x^ value of which was 
determined separately using Monte Carlo simulations. The best 
fit values obtained (x) are a one-dimensional pairwise velocity 
dispersion of o-p = 410±170 km s'^ , and /3 = 0.38±0.17. The 1- 
a and 2-cr error countours are also shown for the fit to the average 
mock catalogue power spectrum (dashed lines) , together with the 
best fit values (+) of ap = 540 km s~\ and /3 = 0.37. 



4 DISCUSSION 

Assuming an EdS universe, we find values for the one- 



dimensional pairwise velocity dispersion of ap 



410 ± 



170 km s~\ and /3 = 0.38 ± 0.17. An open or lambda- 
dominated universe yield similar results (Table nl). x^ con- 
tours for the fit are shown in figure m Although the quoted 
errors assume independence of the two parameters, they are 
clearly related, and a tighter joint constraint would be pos- 
sible. 

It is possible that the assumed P^{k) in the model may 
differ from the true P^(fc) of the Durham/UKST Galaxy 
Survey. Any differences could introduce systematic errors to 
the measurements of ap and /3. To investigate this we re- 
peated the fits described above, simultaneously fitting the 
amplitude of the real-space power spectrum, whilst still as- 
suming the shape of the APM P^{k). The additional pa- 
rameter slightly increased uncertainty in the fitted values of 
ffp and /3, and led to a marginally higher value for the latter 
(/? = 0.51 ± 0.24 in an EdS universe). The fitted amplitude, 
however, was found to be consistent with that of the original 
APM measurement, and we conclude that the assumption 
of the APM P^(fc) anrplitude and shape is adequate for the 
purposes of this paper. However, in a more extensive study, 
using large redshift survey datasets such as 2dF or Sloan, it 
may be necessary to fit the real-space power spectrum self- 
consistently within the model to minimize any systematic 
errors. 

These results are consistent with other estimates from 
the Durham/UKST Galaxy Redshift Survey. Ratcliffe et al. 
(1998a) determined values of /3 = 0.48 and ap = 416 us- 
ing the two-point galaxy correlation function, ^{a, n). Hoyle 
et al. (1999) compared the Durham/UKST redshift-space 
power spectrum to the real-space APM power spectrum in 
order to estimate the redshift distortions, deriving slightly 
higher values of /3 = 0.60 ± 0.35 and ap ^ 450 ± 200 km s"^ 
(Hoyle et al. actually estimated the one-particle velocity dis- 
persion ~ ap/^/2). The /3 estimate favours either a low- 
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